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ABSTRACT 

The  primary  goals  of  this  thesis  included  developing  a  computer-integrated  LDV 
system,  utilizing  the  LDV  system  to  measure  the  mean  and  fluctuating  quantities  of  the 
flow  field  of  an  airfoil-spoiler  system,  and  comparing  the  results  of  the  experiment  with 
previous  hot-wire  anemometry  results. 

The  experiment  was  conducted  in  the  NASA  Ames  1 1  cm  x  25  cm  indraft 
tunnel.  A  two-component  LDV  system  was  developed  and  set  up,  and  measurements 
were  obtained  along  the  upper  surface  of  the  airfoil  (ahead  of  and  behind  the  spoiler) 
as  well  as  in  its  wake,  close  to  the  trailing  edge.  Mean  velocity  profiles  and  turbulent 
stresses  were  calculated  and  comparisons  were  made  with  previous  hot-wire 
experiments. 

It  was  shown  that  the  hot-wire  technique  is  not  as  accurate  as  the  LDV  in  highly 
turbulent  regions  of  the  flow  field.  The  hot-wire  results  are  also  not  reliable  in  the  near 
wake  region  due  to  the  inherent  reverse  flow.  In  other,  less  turbulent,  regions  of  the 
flow  field  where  there  is  no  reverse  flow,  good  agreement  between  previous  hot-wire 
data  and  present  experimental  results  was  found. 
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I.  INTRODUCTION 

A.  SPOILERS 

A  spoiler  is  a  plate  on  the  upper  surface  of  a  wing  which  can  be  deflected  to  spoil 
the  flow  and  to  cause  early  separation.  Consequently,  lift  decreases  and  drag  increases. 
Spoilers  are  used  as  effective  aerodynamic  control  surfaces  in  transport  aircrafts.  They 
can  be  symmetrically  deflected  as  speed  brakes  or  as  lift  dumpers  at  touch  down. 
When  spoilers  are  deflected  asymmetrically,  they  can  provide  substantial  amounts  of 
roll  control  in  certain  conditions  of  flight. 

The  flow  field  of  an  airfoil  with  a  deflected  spoiler  is  complex:  including  flow 
separation,  reattachment  and  vortex  shedding  [Ref.  l:p.  1]. 

The  basic  steady  state  structure  of  the  flow  field  of  an  airfoil  with  a  deflected 
spoiler  has  been  measured  and  hypothesized  by  the  Boeing  Research  Group  [Ref.  2]. 
Further  work  into  how  the  unsteady  wake  relates  to  the  airfoil  geometry  and  to  other 
parameters  was  undertaken  by  Lee,  et  al.  In  their  experiments,  Lee,  et  al.,  used 
hot-wire  anemometers  and  pitot  tubes  to  systematically  measure  the  mean  and 
fluctuating  quantities  on  the  airfoil  surface,  as  well  as  in  the  wake,  in  order  to  better 
understand  the  structure  of  this  flow  field.   [Ref.  1] 

B.  LASER  DOPPLER  VELOCLMETER  SYSTEMS 

Advanced  research  in  experimental  fluid  dynamics  requires  a  familiarity  with 
sophisticated  measurement  techniques.  Optical  methods  and  in  particular,  the  Laser 
Doppler  Velocimeter  (LDV)  are  now  recognized  as  the  most  reliable  means  for 
performing  measurements  in  complex  turbulent  flows  [Ref.  3]. 

Since  its  inception  in  1964,  the  LDV  technique  has  been  revised,  developed,  and 
tested  by  many  researchers  and  its  capabilities  are  reasonably  well-understood.  The 
method  has  the  advantage  of  performing  instantaneous,  non-intrusive  velocity 
measurements  of  small  particles  suspended  in  the  flow.  These  measurements  are 
relatively  independent  of  fluid  properties  and  the  method  does  not  require  frequent 
calibration.  Furthermore,  its  ability  to  measure  the  flow  without  disturbing  it  is 
especially  attractive.  However,  because  of  the  range  of  disciplines  involved  including 
optics,  electronics,  light  scattering,  and  signal  processing,  the  implementation  of  this 
method  can  be  quite  onerous.  Serious  measurement  errors  can  occur  if  the 
measurement  techniques  are  not  applied  correctly. 
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C.       THESIS  GOALS 

Considering  this  previous  work,  the  goals  of  this  thesis  experiment  were  threefold. 
First,  to  develop  and  set  up  a  two-component  LDV  system  with  complete 
computer-experiment  integration  employed.  Second,  to  utilize  this  LDV  system  to 
measure  the  mean  and  fluctuating  velocities  of  an  airfoil  with  a  deflected  spoiler  flap, 
especially  in  the  reverse  flow  region  and  the  near  wake  where  the  fluctuations  are  large. 
Third,  to  compare  the  LDV  data  obtained  with  the  previous  hot-wire  anemometry 
data. 
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II.  EXPERIMENTAL  ASPECTS 

A.  EXPERIMENTAL  SETUP 

The  experiments  were  conducted  in  the  NASA  Ames  1 1  cm  x  25  cm  indraft 
wind  tunnel.  Figure  2.1  shows  the  test  section  of  this  tunnel.  The  test  section  is 
equipped  with  slotted  walls  and  a  number  of  plenum  chambers  above  and  below  the 
centerline.  These  features  were  used  in  previous  experiments  to  give  the  tunnel 
adaptive  walls.  The  adaptive  wall  concept  was  not  employed  in  this  experiment, 
however  the  slotted  walls  did  provide  a  means  of  reducing  wall  interference.  The  wind 
tunnel  also  is  equipped  with  interchangeable  glass  and  plexiglas  test  section  windows. 
For  these  experiments,  the  glass  windows  were  used  due  to  their  higher  refractive 
qualities  which  enhance  the  LDV  signals. 

The  basic  model  used  was  a  Boeing  advanced  transport  research  airfoil  section 
with  a  6-inch  chord  and  maximum  thickness  ratio  of  11.3%  as  shown  in  Figure  2.2. 
This  airfoil  has  a  spoiler  of  15.5%  chord  length  hinged  at  73%  chord.  The  coordinates 
of  the  airfoil  are  given  in  Table  I.  This  airfoil  was  selected  because  of  the  availability 
of  hot-wire  anemometer  and  pitot  tube  experimental  data  for  comparison  with  the 
two-component  LDV  results. 

B.  INSTRUMENTATION 

1.  Pressure  Instrumentation 

For  this  experiment,  pressure  information  was  obtained  by  taking  advantage 
of  a  bank  of  six.  24-port  Scanivalves  tapped  into  the  tunnel  during  a  previous 
experiment.  The  six  pressure  transducers  were  of  the  differential  type  with  a  range  of 
±5  psig.  One  side  of  each  transducer  was  connected  with  the  measured  pressure 
through  plastic  tubing  while  the  other  side  of  each  transducer  was  vented  to 
atmosphere. 

In  order  to  calibrate  the  transducers,  one  port  of  each  Scanivalve  was  tapped 
to  a  calibration  pressure  generated  by  a  small  electric  vacuum  pump.  Another  port  of 
each  Scanivalve  was  tapped  to  the  atmosphere.  The  calibration  pressure  and 
atmospheric  pressure  were  each  measured  by  alternate  methods.  These  two  known 
pressures  and  the  values  given  for  them  by  each  transducer  allowed  the  Scanivalves  to 
be  calibrated. 
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Figure  2.1     NASA  Ames  indraft  wind  tunnel  test  section. 
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Figure  2.2    Boeing  advanced  transport  research  airfoil  section. 

2.    Laser  Doppler  Velocimeter  Instrumentation 

Figure  2.3  shows  the  two-component,  two  color  LDV  arrangement  used  at 
NASA  Ames.  The  laser  is  a  four  Watt  Argon  ion  unit  (Spectra  Physics  Model  64). 
This  laser  is  configured  with  an  end  cavity  multiline  reflector  which  allows  it  to  operate 
in  the  multiline  mode.  Translation  stages  in  three  axes  permit  positioning  of  the  probe 
volume  anywhere  in  the  test  section. 

In  the  multiline  mode,  the  resultant  beam  is  comprised  of  all  wavelengths 
inherent  to  the  Argon  ion  unit.  This  beam  exits  the  laser  in  a  fairly  collimated  state, 
hence  an  external  collimator  is  not  required.  The  beam  first  passes  through  a 
half-wave  plate  (polarization  rotator)  where  it  becomes  horizontally  polarized.  A  set  of 
Brewster  angle  prisms  is  then  used  to  separate  the  primary  beam  into  its  various 
wavelengths.  To  obtain  maximum  light  transmission  through  the  prisms,  the  beam 
polarization  must  be  parallel  to  the  plane  defined  by  the  incident  and  refracted  beams 
[Ref.  4J.  This  was  the  reason  that  the  initially  vertically-polarized  laser  beam  was 
passed  through  a  polarization  rotator  prior  to  entering  the  dispersion  prism.  When  the 
two  prism  method  is  used,  the  prisms  also  turn  the  beams  which,  with  a  little  planning, 
may  be  used  to   advantage.    At   some  distance   from  the   prisms,   the   beams  have 
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TABLE  I 

AIRFOIL  COORDINATES  (INCHES) 

X 

Y 

upper 

Y 

lower 

0.0 

0.0 

0.0 

.0030 

.0222 

-.0108 

.0060 

.0300 

-.0162 

.0150 

.0468 

-.0258 

.0300 

.0672 

-.0348 

.0900 

.1242 

-.0588 

.1200 

.1452 

-.0672 

.1500 

.1626 

-.0750 

.2250 

.2004 

-.0912 

.3000 

.2304 

-.1050 

.4500 

.2790 

-.1296 

.6000 

.3126 

-.1524 

.7500 

.3396 

-.1728 

.9000 

.3594 

-.1920 

1.2000 

.3858 

-.2250 

1.5000 

.3990 

-.2502 

1.S0O0 

.4044 

-.2670 

2.1000 

.4038 

-.274S 

2.4000 

.4002 

-.2742 

2.7000 

.3924 

-.2658 

3.0000 

.3816 

-.2502 

3.3000 

.3678 

-.2298 

3.6000 

.3510 

-.2064 

3.9000 

.3288 

-.1818 

4.2000 

.3018 

-.1560 

4.5000 

.2700 

-.1308 

4.8000 

.2304 

-.1044 

5.1000 

.1842 

-.0792 

5.4000 

.1308 

-.0540 

5.7000 

.0696 

-.0282 

6.0000 

.0024 

-.0024 

sufficient  separation  so  mirrors  can  be  used  to  direct  the  individual  beams  to  the  next 
optical  components.  The  two  stronger  lines,  green  (514.5  nm)  and  blue  (488  nm),  are 
guided  by  means  of  turning  mirrors  to  Bragg  cells.  The  weaker,  unused  lines  are 
"trapped"  (blocked  by  nonreflective  plates)  to  keep  them  from  interfering  with  the 
system.  Prior  to  reaching  its  Bragg  cell,  the  green  beam  is  passed  through  another 
half-wave  plate  in  order  to  vertically  polarize  it. 

A  Bragg  cell  is  an  acousto-optic  device.  Both  the  blue  and  green  beams  are 
passed  through  Bragg  cells  which  split  each  into  two  equal-intensity,  slightly  diverging 
beams  while  inducing  a  40  MHz  modulation  in  one  of  the  beams  in  each  pair.  The 
difference  is  that  the  green  beam  is  split  in  the  vertical  plane  while  the  blue  beam  is 
split  in  the  horizontal  plane  (Figure  2.3). 
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Figure  2.3     Two-component  LDV  arrangement  used  at  NASA  Ames. 
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In  each  case  (blue  and  green),  the  two  diverging  beams  are  then  made  parallel 
by  passing  each  pair  through  its  own  alignment  prism.  From  there  on,  turning  mirrors 
guide  the  beams  through  the  translation  stages.  They  are  finally  focused  by  means  of  a 
condensing  lens  at  the  desired  point  in  the  test  section  of  the  wind  tunnel. 

The  four  beams  must  all  focus  to  a  common  point  in  the  flow  and  be 
approximately  coaxial.  Using  the  steering  mirrors  for  the  appropriate  beam  pair,  the 
focused  beams  are  adjusted  to  visually  overlap.  Following  this,  various  techniques 
exist  for  making  fine  adjustments  to  insure  the  beams  all  cross  at  a  common  point. 

An  off-axis,  forward  scatter  configuration  has  been  selected  for  the  receiving 
optics  used  at  NASA  Ames.  This  gives  flexibility  to  the  system  so  it  can  be  adjusted  to 
different  requirements;  if  maximum  signal  strength  is  required,  angular  deviation  from 
the  beam  axis  of  symmetry  is  minimized.  If  focusing  on  the  smallest  part  of  the  probe 
volume  is  of  interest,  the  angular  deviation  is  made  as  large  as  the  physical  constraints 
will  allow. 

The  first  lens  in  the  receiving  optics  is  the  aperture  through  which  the 
scattered  light  is  collected  (Figure  2.4).  The  second  lens  collimates  the  beam.  A 
collimated  beam  is  required  by  the  polarization  beamsplitter  for  maximum  efficiency. 
The  polarization  beamsplitter  allows  the  vertically  polarized  green  beam  to  pass 
straight  through  while  the  horizontally  polarized  blue  beam  is  reflected  upward.  The 
use  of  polarization  rather  than  frequency  to  separate  the  two  wavelengths  substantially 
increases  the  efficiency  of  the  system.  Each  beam  then  passes  through  a  line  filter 
centered  on  the  respective  wavelength  (blue  or  green)  in  order  to  reduce  both  the 
ambient  light  and  the  cross-talk  due  to  impure  polarization.  Each  line  filter  is  followed 
by  a  lens  which  serves  to  focus  the  light  onto  the  pinhole  aperture  of  the  respective 
photomultiplier.  The  photomultipliers  convert  this  light  energy  into  electronic  signals 
which  are  then  passed  on  to  the  signal  processors. 

Preamplifiers  are  located  in  the  photodetector  housings  to  further  increase  the 
signal  amplitude  and  drive  the  signal  over  relatively  long  cables.  The  preamplifier  is 
kept  close  to  the  detectors  in  order  to  increase  the  signal  level  before  additional  noise 
enters  the  signal  cables.  It  is  always  good  practice  to  amplify  the  signal  as  much  as 
possible  with  the  photomultiplier  tube  before  going  to  the  preamplifier  which  inevitably 
produces  electronic  noise. 

The  signal  processors  used  in  this  experiment  were  of  the  counter  type  and 
produced  by  Macrodyne.    These  processors  have  good  resolution  and  superior  noise 
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Figure  2.4    Computer-experiment  integration. 
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rejection  circuitry.  The  M aerodyne  unit  is  a  high-speed  counter  with  a  1  GHz  clock 
designed  to  operate  in  a  variety  of  experimental  situations.  It  has  been  applied  under 
single  particle,  multi-particle  and  multi-velocity  conditions.  In  addition  to  the  basic  5/8 
(or  10/16)  signal  validation  criteria  and  oversize  particle  rejection,  a  sophisticated 
multi-level  sequence  validation  has  been  incorporated.  This  allows  the  experimenter  to 
select  the  criteria  which  will  best  yield  data  for  a  given  signal-to-noise  ratio,  signal 
level,  seeding  level,  Doppler  frequency,  and  fluctuation,  indicative  of  the  test  situation. 

When  the  particle  seeding  rate  is  low,  as  was  the  case  in  this  experiment,  the 
LDV  signals  become  intermittent  bursts.  Under  these  conditions,  the  counter 
processor  is  the  best  choice  as  a  means  for  obtaining  the  Doppler  difference  frequency. 
If  the  seeding  is  high,  the  signals  will  be  almost  continuous  with  random  amplitude  and 
phase  fluctuations.  The  counter  will  still  function  but  cannot  be  expected  to  produce 
accurate  results  to  better  than  ±  1/2  cycle. 

The  counter  processor  determines  the  Doppler  difference  frequency  in  the  time 
domain.  Timing  over  a  selected  number  of  cycles  in  the  burst  signal  is  accomplished 
with  a  crystal  controlled  "clock"  oscillator.  The  frequency  of  the  oscillator  must  be 
much  higher  than  the  Doppler  difference  frequency  to  achieve  good  resolution.  The 
clock  count  is  initiated  when  the  high-pass  filtered  signal  exceeds  a  set  threshold  level. 
A  preselected  number  of  cycles  (8  or  16)  are  processed  and  used  to  stop  the  counter. 
Greater  accuracy  may  be  obtained  by  processing  more  cycles.  Some  processors,  (e.g., 
Aerometrics)  count  all  of  the  cycles  in  the  burst  and  average  the  result  over  the  number 
of  cycles  counted.  This  capability  is  especially  important  when  the  seed  particles  are 
polydisperse  and  produce  a  range  of  signal  amplitudes.  In  such  cases,  there  is  the 
possibility  of  completing  the  preset  number  of  cycles  on  the  early  part  of  the  signal 
which  has  a  lower  signal-to-noise  ratio  and  slew  rate.  There  is  also  the  possibility  of 
reading  signals  from  larger  particles  more  than  once. 

Although  the  importance  of  attaining  good  quality  signals  cannot  be 
overemphasized,  in  many  practical  cases  the  signals  will  be  noisy.  Therefore,  each 
burst  must  pass  certain  criteria  before  it  is  accepted.  The  so-called  periodicity  check  is 
used  to  reject:  1)  signals  with  insufficient  cycles,  2)  signals  produced  by  two  or  more 
particles  in  the  measuring  volume  at  a  time,  and  3)  noise.  Two  counters  are  used  to 
make  period  measurements  over  a  different  number  of  cycles,  usually  5  and  8.  The 
ratio  of  the  two  measurements  is  then  compared  to  five-eighths  and  if  they  agree  to  a 
specified  percentage,  the  signal  is  validated.    The  problem  with  this  technique  is  that 
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the  entire  signal  may  be  rejected  because  of  a  noisy  cycle  early  in  the  burst.    Howe   . 
when  the  method  is  coupled  with  the  three-level  validation  circuitry,  it  can  be  effective 
in  rejecting  spurious  noise. 

The  three-level  validation  requires  that  a  signal,  Figure  2.5.  exceed  a  positive 
threshold  level,  pass  through  zero,  and  exceed  a  negative  threshold  lei ..  bel  :re  the  next 
zero  crossing  is  counted  as  a  cycle.  Such  logic  has  proven  to  be  very  effective  in 
rejecting  noise.  If  the  conditions  are  not  met  on  any  cycle,  the  counter  resets  and 
starts  again  almost  immediately.  Thus,  noisy  cycles  on  the  start  of  a  burs:  need  not 
.    use  the  loss  of  an  otherwise  good  signal. 

These  amplified  electronic  signals  are  sent  to  the  signal  processor  and 
ultimately  to  the  data  processor  for  conversion  into  velocities. 

Ideally,  an  LDV  needs  no  calibration  since  it  provides  a  direct  measure  of  the 
panicle  velocity.  The  fringe  spacing.  6,  can  be  determined  using  the  following  formula 
where  /.  is  the  laser  wavelength  and  9  is  the  angle  of  the  intersecting  beams: 

5  =  V(2  *  sir..  0  2 

Knowing  the  fringe  spacing  and  the  time  required  for  a  given  number  of  cr:  gs 
particle  velocity  can  be  determined.  In  practice  this  method  can  be  inaccurate  due  tc 
misalignment  of  the  optics.  The  beam  intersection  angle,  0.  is  also  very"  difficult  tc 
compute  accurately.  Therefore,  it  is  wise  to  test  the  entire  system  by  measuring  the 
light  scattered  from  an  object  moving  at  a  velocity  that  is  known  a  priori.  A  rota:::-.. 
disk  with  a  mark  at  a  known  radius  is  usually  the  most  convenient  means  of  doing  this. 
However,  another  method  (the  method  used  in  this  experiment)  is  to  measure  the  free 
stream  velocity  using  the  LDV  and  compare  this  value  with  that  computed  from  the 
Scaniva/.  e  data. 

3.   Flow  Seeding  Requirements 

An  inconvenient  characteristic  of  the  laser  Doppier  velocimeter  techi  .e  is 
-  need  for  seed  particles.  The  LDV  actually  measures  the  velocity  of  small  particles 
that  are  assumed  to  move  with  the  flow.  If  these  particles  do  not  have  the  appropriate 
size  and  concentration  distribution,  the  experiment  results  may  be  adversely  affected 
Particle  sizes  on  the  order  of  1  micron  are  often  quoted  as  adequate  for  mos:  gas  (lows 
and  can  respond  to  turbulence  frequencies  exceeding  one  kHz.    Smaller  particles  ma] 
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Figure  2.5    Three-level  detection  scheme. 
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be  below  the  detectable  limit  of  the  instrument  but  will  contribute  to  the  background 
noise  while  larger  particles  that  are  easily  detectable  may  not  adequately  track  the  flow. 
It  was  determined  in  this  experimental  situation  that  seeding  of  the  flow  was 
not  absolutely  necessary.  Several  runs  were  made  with  no  seeding  supplied.  The  data 
obtained  in  these  runs  was  accurate  as  dust  particles  acted  as  the  flow-following  agent. 
However,  it  was  determined  that  by  using  a  flow  seeding  system  the  data  rate  was 
increased  by  a  factor  of  four  to  five.  The  system  used  consisted  of  a  paint  sprayer 
spraying  latex  particles  (on  the  order  of  3  microns)  suspended  in  alcohol. 

C.  COMPUTER-EXPERIMENT  INTEGRATION 

The  experiment  was  fully  integrated  with  a  Data  General  NOVA  minicomputer 
for  on-line  data  acquisition  and  analysis,  Figure  2.4.  A  main  program 
("SPOILER.FR")  was  written  to.  perform  all  tasks  and  record  all  data  required  during 
the  experiment  so  that  individual  tunnel  runs  proceeded  virtually  'hands  ofT.  The  main 
program  called  a  host  of  other  programs  during  each  run  in  order  to  accomplish 
various  tasks.   These  other  routines  are  listed  here: 

(1)  "SCANT. FR"  -  used  to  obtain  pressure  data  from  the  Scanivalves  and 
compute  free  stream  velocity  and  free  stream  Mach  number.  "SCANT. FR" 
itself  calls  three  other  routines  in  order  to  determine  ambient  temperature, 
calibration  pressure  and  barometric  pressure: 

(a)  "TEMP.FR"  -  which  vields  ambient  temperature  from  a  thermocouple 
via  an  analog-to-digital  converter.  The  thermocouple  is  located 
upstream  of  the  test  section. 

(b)  "CALPRES.FR"  -  which  calculates  the  output  pressure  from  the 
calibration  vacuum  pump.  The  pressure  value  is  read  by  this 
program  using  an  analoe-to-digital  converter.  This  value  is  then  used 
in  the  calibration  of  the  scanivalves. 

(c)  "NPARO.FR"  -  which  yields  the  value  for  current  atmospheric 
pressure  from  a  Paroscientific  pressure  sensor  via  a  digital  I/O. 

(2)  "STDTM.FR"  -  which  gives  the  current  time  and  date  in  order  that  the  data 
files  being  created  by  the  main  program  can  be  labeled. 

(3)  "TRELAY.FR"  -  to  command  the  taskmaster  to  move  the  traverse  mechanism 
to  the  next  data  collection  location. 

(4)  "PLOTP.FR"  -  to  provide  on-line  data  analysis  by  allowing  any  of  five  curves 
to  be  plotted  immediately  after  a  given  run. 

D.  EXPERIMENTAL  CONDITIONS 

The  free  stream  velocity  was  maintained  at  a  nominal  value  of  70  m/sec, 
corresponding  to  a  Reynolds  number  of  7.4  x  105  based  on  the  airfoil  chord. 
Experiments  were  carried  out  with  the  spoiler  deflected  30  degrees  and  the  angle  of 
attack  at  0  and  5  degrees.    Locations  where  data  was  obtained  are  shown  in  Figure  2.6. 
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Figure  2.6    Data  collection  locations. 
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In  order  to  determine  the  velocity  profiles,  data  was  acquired  at  varying  steps 
based  on  the  velocity  gradient.  Fine  step  sizes  were  taken  near  the  airfoil  for  boundary 
layer  profiles  and  near  free  shear  layers  in  the  wake  measurements.  Typically,  15  data 
points  were  taken  in  the  boundary  layer  and  20  points  in  the  wake.  The  4-beam  LDV 
arrangement  precluded  obtaining  data  closer  than  0.1  inch  from  the  airfoil  surface, 
however. 

The  boundary  layer  traverse  was  along  the  vertical  axis  of  the  test  section,  so  the 
traverse  direction  is  not  normal  to  the  airfoil  surface  as  are  usual  boundary  layer 
profiles.  Angle  of  attack  and  model  surface  curvature  make  the  measurement  axis 
incline  to  the  model. 
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III.  DATA  ACQUISITION 

A.       FREE  STREAM  VELOCITY 

At  the  beginning  of  every  run,  the  main  experiment-driving  program 
("SPOILER. FR")  calls  for  the  Scanivalve  pressures,  ambient  temperature,  calibration 
pressure  and  atmospheric  pressure  to  be  read.  Assuming  isentropic  flow  in  the  test 
section  upstream  of  the  airfoil  and  away  from  any  shock  waves,  the  free  stream  Mach 
number,  M,  can  be  calculated  from  the  stagnation  pressure,  P  ,  and  the  static 
pressure,  P,   along  the  wall: 

M  =  {((P0/P)^-l)/Y-i)(2/(y-l))}l/2 
For  y—  1.4: 

M  =  {((P0/P)^2/7)-l)x5}1/2 
Then,  velocity  is  given  by: 

V  =  Ma 


where 


a  =  vTRT" 
T  is  obtained  from  the  total  temperature,  TQ,  using  the  following  formula: 

T  =  T0/{1  +  ((y-l)/2)xM2} 
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TQ  is  measured  by  a  thermocouple  at  the  wind  tunnel  inlet. 

The  above  method  is  used  to  obtain  free  stream  Mach  number  and  free  stream 
velocity  for  every  run.  The  value  for  free  stream  velocity  is  used  in  the  data  reduction 
process  to  normalize  the  airfoil  velocities  obtained  from  the  LDV. 

B.        AIRFOIL  VELOCITIES 

Processing  of  the  Laser  Doppler  Velocimeter  data  obtained  during  this 
experiment  was  carried  out  using  the  same  Data  General  Nova  digital  computer  and 
driving  program  used  for  all  other  experiment  functions.  "SPOILER. FR",  the  main 
program  for  this  experiment,  calls  a  digital  I/O  driver  program  ("DI024")  to  read 
words  as  fast  as  possible  from  two  cards  (these  cards  receive  the  digital  output  from 
the  two  VI aerodyne  channels,  one  for  each  LDV  component).  As  soon  as  1000 
samples  (per  component)  have  been  taken  by  the  Macrodyne  processor  and  read  by 
DI024,  the  binary  output  from  each  sample  is  converted  to  decimal  using  a  decoding 
subroutine  ("DECODE l.FR"). 

As  mentioned  earlier,  the  Macrodyne  signal  processor  used  in  this  experiment 
utilizes  a  1  GHz  counter.  The  number  of  "counts"  required  for  the  passage  of  8  cycles 
is  sent  to  the  digital  I/O  in  the  form  of  a  10-bit  mantissa  and  4-bit  exponent,  namely 

D9,  D8,...D0  x  2exP 

For  example,  if  400  "counts"  were  required  for  a  particular  Doppler  burst,  the 
output  would  be 

0000110010      0011 
Mantissa         EXP 


which,  upon  decoding,  reduces  to 

50  x  23  =  400 
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Since  the  clock  counts  are  produced  at  4  GHz  (the  1  GHz  clock  is  divided  into  4 
phases),  each  clock  count  is  0.25  nanoseconds.   Thus, 


exn           °-25) 
t.  =  Dg,  D„,...Dn  *  2exP  x  ns 

S  9        8  0  (80) 


is  the  period  of  the  signal  when  8  cycles  are  measured. 
Converting  this  to  signal  frequency,  fs,  yields 


1  25"exP  x  1000 

f _  =  —  =  MHz 

t  Mantissa 


Since   the   signal   frequency,   f§,    is   composed    of  the   Doppler   frequency,    the   shift 
frequency  and  the  mixing  frequency,  the  Doppler  frequency  can  be  obtained  from 


fD  "  fs       fshift  +  fmix 


Once  the  Doppler  frequency  of  a  particle  is  known,  its  velocity  can  be  determined 
from 


velocity  =  const  x  fi 


D 


where  "const"  is  determined  based  on  the  fringe  spacing  or  direct  calibration. 

With  the  1000  samples  from  each  component  converted  to  decimal,  the  data  can 
be  reduced  to  average  and  rms  velocities  using  the  formulas  found  in  Table  II.  This 
data  (both  raw  and  reduced)  is  then  stored  for  future  analysis. 

While  most  data  analysis  is  performed  after  completion  of  the  experiment,  several 
on-line  analysis  routines  are  available  for  use  as  real  time  performance  checks.  For 
example,  "HISTOGRAM. FR"  is  a  program  available  to  first  read  a  given  number  of 
samples  from  one  of  the  Macrodyne  channels  and  plot  frequency  versus  the  number  of 
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samples  at  each  frequency  in  histogram  fashion.  With  this  program  it  is  easy  to 
determine  whether  or  not  the  data  is  conforming  to  the  Normal  Distribution  (to  which 
ideal  data  should  conform).  Another  on-line  analysis  tool  which  is  built  into  the  main 
program  "SPOILER. FR"  is  the  subprogram  "PLOTP.FR".  This  routine  gives  the 
experimenter  the  option  to  have  any  of  five  velocity  profiles  (i.e.,  Uave/Uooi 
Vave/V00,  etc.)  plotted  on  the  screen  immediately  after  completing  an  experimental 
run.  This  subprogram  was  used  extensively  to  immediately  inspect  the  quality  of  each 
run. 


TABLE  II 

VELOCITY  PROFILE  AND  TURBULENT  STRESS  EQUATIONS 
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C.   REYNOLDS  SHEAR  STRESSES  AND  TURBULENT  STRESSES 

Table  II  shows  the  equations  used  to  determine  mean  velocities,  turbulent 
stresses  and  Reynolds  shear  stress  throughout  this  experiment.  The  value  for  N  was 
generally  1000  since  this  was  the  number  of  velocity  samples  taken  at  each 
measurement  location  (for  each  component). 
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IV.  RESULTS  AND  DISCUSSION 

A.  INTRODUCTION 

One  of  the  purposes  of  this  experiment  was  to  demonstrate  the  benefits  of  the 
two-component  Laser  Doppler  Velocimeter  by  employing  this  type  of  system  to 
measure  the  flow  field  of  an  airfoil  with  a  deflected  spoiler.  In  order  to  present  the 
results  of  these  measurements  effectively,  a  graphics  program  called  "PLOT-3D"  was 
used  to  plot  the  data  obtained.  "PLOT-3D"  is  available  on  the  VAX  system  at  the 
NASA  Ames  Research  Center. 

A  data  interpolation  routine  written  by  Chyang  Sheng  Lee  of  Stanford  University 
was  also  used  to  interpolate  the  measured  data  into  a  grid  system  so  that  it  could  be 
plotted.  All  of  the  contour  plots  presented  use  the  interpolation  routine  to  interpolate 
in  both  the  x  and  y  directions.  However,  some  of  the  velocity  vector  plots  use 
interpolation  in  the  y  direction  only.  This  was  done  because,  in  some  cases,  linear 
interpolation  of  the  data  gives  a  false  impression  of  what  the  flow  field  is  actually 
doing. 

In  all  cases,  the  data  presented  in  this  section  is  for  a  =  0  :  and  5  =  30°. 

B.  FLOW  FIELD  ANALYSIS 

Velocity  measurements  were  taken  at  a  variety  of  points  along  the  upper  surface 
ahead  of  and  behind  the  spoiler.  Figure  4.1  shows  a  plot  of  the  data  obtained  over  the 
rear  portion  of  the  airfoil.  Note  the  low-velocity,  recirculating  flow  region  dire.:/. 
behind  the  spoiler. 

While  the  data  for  Figure  4.1  was  interpolated  in  the  y  direction  only,  the  data 
for  Figure  4.2  was  interpolated  in  both  the  x  and  y  directions.  It  can  be  seen  in  Figure 
-.2  that  interpolation  in  the  x  direction  yields  some  fairly  unrealistic  velocity  vectors 
due  to  the  relatively  low  number  of  measurement  locations  as  compared  to  the  number 
of  grid  points.  For  example,  those  directly  ahead  of  the  spoiler.  On  the  other  hand, 
complete  interpolation  (in  both  the  x  and  y  directions)  does  serve  to  paint  a  good 
overall  picture  of  what  the  flow  field  is  doing  based  on  a  finite  number  of  velocity 
measurer/.er.is. 
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Figure  4.1     Flow  field  measured  over  rear  portion  of  airfoil 
(interpolation  in  y  direction  only). 
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Figure  4  2     Flow  Held  measured  over  rear  portion  of  airfoil 
(interpolation  in  both  directions). 
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C.       WAKE 

The  wake  behind  the  spoiler  was  measured  with  the  two-component  LDV  system 
at  numerous  locations  to  determine  the  mean  and  rms  profiles.  The  results  from  the 
wake  are  presented  in  this  section. 

1.  Wake  Profiles 

The  mean  longitudinal  and  transverse  wake  profiles  measured  by  the  LDV  for 
a  =  0°  and  6  =  30°  are  presented  in  Figures  4.3  and  4.4.  The  symbols  are  the 
measured  data  and  the  lines  represent  the  best  fit  curves  of  the  data.  As  expected,  the 
width  of  the  wake,  the  velocity  defect  and  transverse  velocity  all  decrease  in  the 
downstream  direction. 

The  Reynolds  shear  stress,  longitudinal  and  transverse  turbulent  stresses 
measured  by  the  LDV  are  presented  in  Figures  4.5  to  4.7.  Obvious  in  all  three  of  these 
figures  is  the  fact  that  the  fluctuating  region  appears  to  grow  from  the  x/C=1.0 
(trailing  edge)  location  to  x/C=  1.25  then  decrease  as  it  approaches  x/C=  1.5.  The  "S" 
shape  profile  in  the  Reynolds  shear  stress  is  typical  for  a  two-dimensional  wake 
[Ref  l:p.  29].  The  longitudinal  and  transverse  turbulent  stresses  show  two  peaks  at  the 
locations  where  the  two  free  shear  layers  (corresponding  to  the  airfoil  and  the  spoiler) 
oscillate  in  the  near  wake. 

Figures  4.8  and  4.9  show  the  velocity  vectors  in  the  wake.  A  scaled  profile  of 
the  airfoil  is  superimposed  in  each  figure  for  reference.  The  data  in  Figure  4.8  was 
interpolated  in  the  y  direction  only.  Figure  4.9  was  interpolated  in  both  directions  and 
depicts  the  mean  wake  velocity  profile  quite  well. 

2.  Contour  Plots 

Figures  4.10  to  4.16  depict  "contour"  plots  of  various  quantities  in  the  airfoil 
wake.  Figure  4.10  shows  regions  of  equivalent  U-magnitude  (U  is  referenced  to  free 
stream  velocity).  Note  that  there  are  regions  where  the  ratio  U/Uqo  *s  greater  than 
unity.  This  is  due  to  tunnel  blockage  effects  which  cause  the  flow  to  accelerate  more 
than  normal  around  the  airfoil. 

Figure  4.11  depicts  the  regions  of  equivalent  V-magnitude  (also  referenced  to 
free  stream  velocity).  A  minus  sign  on  a  contour  line  in  Figure  4. 1 1  implies  that  the  V 
component  of  velocity  is  negative  (downward)  in  that  region. 

In  Figure  4.12,  the  total  velocity  magnitudes  in  the  wake  are  contoured. 
Again,  contour  lines  with  values  greater  than  one  arise  due  to  blockage  effects. 
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Figure  4.3     Longitudinal  velocity  profiles  in  the  wake. 
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Figure  4.4    Transverse  velocity  profiles  in  the  wake. 
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Figure  4.5     Reynolds  shear  stress  in  the  wake. 
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Figure  4.6     Longitudinal  turbulent  stress  in  the  wake. 
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Figure  4.7    Transverse  turbulent  stress  in  the  wake. 
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Figure  4.8    Velocity  vectors  in  the  wake 
(interpolation  in  y  direction  only). 
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Figure  4.9     Velocity  vectors  in  the  wake 
(interpolation  in'  both  directions). 
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Figure  4. 10    Contour  plot  of  U-magnitude  in  the  wake. 
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Figure  4. 1 1     Contour  plot  of  V-magnitude  in  the  wake. 


43 


CO 


Figure  4.12    Contour  plot  of  total  velocity  magnitudes 
in  the  wake. 
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Figure  4.13     Contour  plot  of  flow  angles  in  the  wake. 
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Figure  4.14    Contour  plot  oflonsitudinal  turbulent  stresses 

in  the  waRe. 
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Figure  4.15     Contour  plot  of  transverse  turbulent  stresses 

in  the  wake. 
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Figure  4.16    Contour  plot  of  Reynolds  shear  stresses 
in  the  wake. 
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Figure  4.13  shows  how  rapidly  the  flow  angles  change  in  the  near  wake  (the 
closer  the  contour  lines,  the  steeper  the  gradient). 

Longitudinal  and  transverse  turbulence  and  Reynolds  shear  stress  are 
contoured  in  Figures  4.14  to  4.16.  Figure  4.14  clearly  shows  the  two  regions  where  the 
turbulence  intensity  is  the  highest.  These  two  regions  are  probably  caused  by 
counter-rotating  vortices  originating  at  the  spoiler  trailing  edge  (clockwise  vortex)  and 
the  airfoil  trailing  edge  (counterclockwise  vortex). 

D.       COMPARISON  WITH  HOT-WIRE  ANEMOMETER  DATA 

One  purpose  of  this  thesis  was  to  compare  the  data  obtained  in  this  experiment 
with  the  hot-wire  data  obtained  under  similar  conditions  by  Lee  [Ref.  5].  Figures  4.17 
to  4.21  are  LDV/hot-wire  comparisons  of  longitudinal  and  transverse  velocity, 
longitudinal  and  transverse  turbulent  stress,  and  Reynolds  shear  stress  at  three 
locations  in  the  wake. 

It  is  important  to  mention  that,  in  the  previous  experiment,  no  hot-wire  data  was 
obtained  in  the  region  of  reverse  flow  between  x/C=1.0  and  x/C=1.25.  Hot-wire 
anemometry  is  virtually  unusable  in  this  region  due  to  the  fact  that  hot-wires  rectify  the 
signals,  which  include  the  instantaneous  forward  velocity  as  well  as  the  reverse  velocity. 
Consequently,  velocities  measured  by  hot-wires  in  reverse  flow  regions  appear  higher 
than  they  actually  are. 

Except  for  minor  grid  alignment  errors  in  the  y  direction,  most  of  the  curves  are 
in  very  good  agreement.  As  expected  however,  the  longitudinal  turbulence  stress 
curves  at  x/C=1.25  (Figure  4.19)  and  the  transverse  turbulence  stress  curves  at 
x/C=1.25  (Figure  4.20)  show  considerable  disagreement.  In  the  longitudinal 
turbulence  case,  the  hot-wire  data  shows  more  turbulence  in  the  middle  of  the  wake 
while  in  the  transverse  turbulence  case,  the  hot-wire  data  shows  less  turbulence  in  the 
middle  of  the  wake.  This  is  again  due  to  the  fact  that  hot-wire  measurements  give  a 
misleading  velocity  profile  close  to  the  trailing  edge  because  of  the  reverse  flow  in  that 
region  [Ref.  l:p.  70]. 

The  hot-wire  anemometer  is  quite  accurate  in  regions  where  there  is  no  reverse 
flow.  However,  in  order  to  obtain  reliable  data  on  the  mean  velocities  and  turbulent 
stresses  in  the  reverse  flow  regions  and  the  near  wake  of  this  airfoil-spoiler  system, 
Laser  Doppler  Velocimetry  yields  better  results. 
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Figure  4.17    Comparison  of  LDV  and  hot-wire  data  for 
longitudinal  velocity. 
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Figure  4.18     Comparison  of  LDV  and  hot-wire  data  for 
transverse  velocity. 
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Figure  4.19     Comparison  of  LDV  and  hot-wire  data  for 
longitudinal  turbulence  stress. 
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Figure  4.20    Comparison  of  LDV  and  hot-wire  data  for 
transverse  turbulence  stress. 
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Figure  4.21     Comparison  of  LDV  and  hot-wire  data  for 
Reynolds  shear  stress. 
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E.       ERROR  ANALYSIS 

In  this  experiment,  1000  velocity  samples  were  taken  by  the  LDV  in  each 
direction  (U  and  V)  at  every  measurement  location.  The  purpose  in  taking  so  many 
samples  was  to  provide  enough  values  in  the  fluctuating  regions  to  ensure  that  an 
accurate  statistical  mean  could  be  computed.  Consequently,  in  the  data  reduction 
process  all  1000  values  were  averaged  to  obtain  values  for  Uave  and  V  at  each 
location.  The  equations  in  Table  II  were  then  used  to  obtain  Reynolds  shear  stress, 
longitudinal  and  transverse  turbulent  stress  at  these  locations  (still  using  all  1000 
samples). 

To  determine  its  effect,  a  different  data  processing  method  was  tested  on  the  raw- 
data  from  several  runs.  In  this  procedure,  only  the  samples  which  fell  within  one 
standard  deviation  of  the  mean  were  considered  (approximately  67%  of  the  1000 
samples)  in  computing  L*ave,  V  ,  L'rms  and  Vrms.  Again  the  equations  in  Table  II 
were  employed,  just  as  before.  This  new  process  had  little  effect  on  Lave  but  V 
changed  considerably  (from  -1.3  to  -11.0  at  one  location  tested).  However,  the  V 
values  did  not  always  decrease  (one  location  changed  from  8.7  to  16.9).  Therefore,  the 
changes  were  not  due  to  the  removal  of  a  bias  error.  The  longitudinal  and  transverse 
turbulence  stress  values  decreased  in  all  cases.  This  was  expected  since  the  biggest 
contributors  to  these  values  are  the  velocity  terms  which  differ  greatly  from  the  mean. 

This  new  data  processing  procedure  did  not  improve  the  data  quality  by 
removing  any  biases  (based  on  the  fact  that  some  values  of  L'ave  and  Vave  increased 
while  others  decreased  on  using  the  new  process).  As  a  matter  of  fact,  this  process  had 
a  negative  effect  on  the  data  quality  in  that  it  neglects  velocity  extremes  which,  in  turn, 
lowers  the  apparent  turbulence  intensity  of  the  flow  field.  For  this  reason,  the  original 
data  processing  procedure  (consider  all  1000  samples)  was  used  in  this  experiment. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

A  two-component  LDV  system  was  set  up  for  the  purpose  of  measuring  mean 
and  fluctuating  quantities  of  an  airfoil  with  a  deflected  spoiler.  The  results  of  this 
experiment  have  shown  that  this  system  provides  an  efficient  and  effective  means  of 
flow  field  measurement,  including  regions  of  reverse  flow  and  large  velocity 
fluctuations.  Comparison  with  hot-wire  data  demonstrates  the  LDV's  superiority  over 
the  hot-wire  in  measuring  in  the  near  wake  region. 

The  results  also  clearly  show  two  recirculating  regions  in  the  wake  of  the  airfoil. 
One  vortex  originates  at  the  trailing  edge  of  the  airfoil  itself  while  the  other  originates 
at  the  trailing  edge  of  the  spoiler.  The  turbulence  region  is  concentrated  in  the  near 
wake  and  dissipates  rapidly  in  the  downstream  direction. 

B.  RECOMMENDATIONS 

In  order  to  establish  a  more  complete  data  base,  it  is  recommended  that  more 
LDV  measurements  be  obtained  at  other  than  a  =  0°  and  6  =  30°.  These  results  are 
important  to  the  modeling  of  this  flow  field.  Also,  to  remove  the  possibility  of  tunnel 
blockage  errors,  it  is  necessary  to  perform  further  studies  in  a  large  scale  wind  tunnel. 

When  using  an  LDV  to  measure  a  turbulent  flow  field  such  as  in  this  experiment, 
a  large  number  of  samples  should  be  taken  at  each  location.  Furthermore,  unless 
attempting  to  remove  a  velocity  bias  all  of  the  samples  taken  should  be  considered  in 
computing  Uave,  Vave,  etc. 
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APPENDIX 
COMPUTER  PROGRAMS 


C   PURPOSE  -  (1)  TO  CONVERT  A  24  HOUR  CLOCK  TO  A  12  HOUR  CLOCK. 
C  -  (2)  TO  ADD  A.M.  OR  P.M.  TO  THE  "ADJUSTED"  TIME. 

C 
C 

SUBROUTINE  STDTM  ( I  DATE. ITIME. I VAL) 

D I  THINS  ION  1DATEC3) 

DIMENSION  ITIHEO) 

CALL  DATE  ( 1  DATE , I  ERROR ) 

CALL  TIME  ( IT I  ME . I  ERROR ) 

IF  (ITIMEU  )  .GE.13)  GO  TO  30 

IF  (ITIMEU  )  .LQ.  12)  GOTO  10 

CO  TO  20 
10    IF  (ITIME(2) .GT.0O)  CO  TO  30 

IF  UTIME(3)  .GT.0O)  GO  TO  30 
20    I VAL  =  "AM" 

GO  TO  40 
30    I VAL  =  "PM" 

IF  (ITIMEU  )  .LT. 13)  GO  TO  40 

ITIMEU)  =  ITIMEU)  -  12 

CO  TO  50 
40    IF  (ITIMEU)  .GE. 01  )  GO  TO  50 

ITIMEU)  =  ITIMEU)  +  12 
50    CONTINUE 

RETURN 

END 
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G  ****************************** 
C  **********************-.':*-•'••***** 
c  **  :,:* 

c  **  SP01LER.FR  ** 
c   **  *# 

G  **  MIKE  FOREMAN  ** 
C  **  ** 

G  *****************************:;: 

G  ****«:******#*********#******** 

G 

C  LASER  DATA  ACQUISITION  PROGRAM 

C 

G  RLDll  SPOILER  STDTM  D1024  &FLIB9 


C 


DIMENSION  IFILE(4),  IDT<10),  VAR(IO) 

DIMENSION  ITIIIEO),  IDATEO) 

DIMENSION  IRAYl(lOOO),  1RAY2U000) 

DIMENSION  IDINK16),  1BIN2U6),  VEL1  (  1000  )  ,VEL2(  1  000  ) 

DIMENSION  V(5,20)  ,X(3«./)  ,Y<30) 


G 

C  INITIALIZE 


t. 


DO    10    1=1,3 
IDT( I )=0 
VAR( I )=10. 
10    CONTINUE 

ACCEPT    "RUN    NUMBER    (ZERO  TO  QUIT):     ",IRUN 
IF( IHUN.LE.O)    GO   TO    999 
15    ACCEPT    " NUMBER   OF    SAMPLES    PER   COMPONENT:     " , NUW 
IF(NUM.LE. 1000. AND. NUM. CE.O)    CO   TO    20 


TYPE  "NUMBER  OF  SAMPLES  MUST 

BE  BETWEEN 

CO  TO  15 

ACCEPT  "MIXING  FREQUENCY  #1: 

" , FREfl 1 

ACCEPT  "MIXING  FREQUENCY  *2: 

" , FREQ2 

1    AND    1000" 
CO   TO    15 
20 

C0N1 =16.65 

CON2=16.65 

VAIK  I  )=C0N1 

VAIK2)=C0N2 
C 

C  SUBDEVICE  ADDRESSES 
G 

ISUBl=46K 

ISUB2=47K 
G 

G  TIME  PROGRAM 
G 

CALL  STDTM ( I  DATE , I T I ME , I VAL ) 
G 

DO  30  1=1, 3 

II)T(  I  )  =  1DATE(  I  ) 

1DT( I+3)=1TIME(I ) 
30  CONTINUE 

IDT(7)=IVAL 

IDT(3)=IRUN 

VAR(9)=FREQ1 

VAIK  10)=FREQ2 
C 

35  TYPE  "WILL  THIS  BE  A  :   (1)  WAKE  PROFILE" 

ACCEPT  "  (2)  BOUNDARY  LAYER  PROFILE  " . I  PROF 

G 

IFdPROF.EQ.l  .OR.  I  PROF.  EQ.  2)    GO   TO   40 

TYPE  "MUST  ENTER  1  OR  2" 
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CO  TO  35 

40  IF( I  PROF. EQ. 2)  CO  TO  41 

ACCEPT  "X-LOCATION  FORTIUS  RUN  (XT/C)  :  ",VAR(3) 
CO  TO  42 

41  ACCEPT  "X-LOCATION  FOR  THIS  RUN  (X'C):  ",VAR(3) 

42  ACCEPT  "NUMBER  OF  Y-LOCATIONS:  "  .  NflAX 
IF( I  PROF. EQ. 2)  CO  TO  45 

ACCEPT "STANDARD  WAKE  PROFILE  Y-LOCATIONS?  (l=YES,0=NO)  ",LL 

IF(LL.EQ. 1 )  GO  TO  44 

TYPE" INPUT  Y-LOCATIONS  IN  INCHES  (LOWEST  POINT  FIRST):" 

DO  43  1=1, NNAX 

WRITE  (10.49)  I 

ACCEPT"  ",Y(1) 

43  CONTINUE 
CO  TO  50 

44  Y( 1 )=-l  .0 
Y(2)=-0.8 
Y(3)=-0.6 
Y(4)=-0.46 
Y(5)=-0.34 
Y(6)=-0.2 
Y(7)=-0.1 
Y(8)=0.O 
Y(9)=0.1 
Y( 10)  =  . 2 
Y(ll  )=.3 
Y( 12)= .4 
Y( 13)  =  . 5 
Y( 14)= .6 
Y(15)=.7 
Y( 16)=  .0 
Y(17)=.9 
Y( 10) = 1 .0 
Y( 19)=1 .2 
Y(20)  =  l  .4 
CO  TO  50 

45  TYPE' STANDARD  BOUNDARY  LAYER  PROFILE  Y-LOCATIONS?  ■ 
ACCEPT"  (1=YES,  0=N0):  " , LK 
IFU.K.EQ.  1  )CO  TO  47 

TYPE" INPUT  Y-LOCATIONS  IN  INCHES  (LOWEST  POINT  FIRST):" 
DO  46  1=1 .NMAX 
WRITE (10.49)  I 


ACCEPT"  ",Y(I) 

46  CONTINUE 

CO  TO  50 

47  Yd  )=.  1 

Y(2)= .2 

Y(3)=.3 

Y(4)= .4 

Y(5)=.5 

Y(6)= .6 

Y(7)= .7 

Y(8)  =  .8 

Y(9)= .9 

Y( 10) = 1 .0 

Y( 1  1  )  =  1  .  1 

Y( 12)  =  1  .2 

Y( 13  )  =  1  .3 

Y(14)  =  I  .4 

Y( 15)  =  1  .5 

49  F0RM\T(6X. *Y( " , 12 

"  )  :  "  ) 

50  WRITE( 10. 125) 

[25  FORMAT (  "  "  ,  "BEGIN 

(SCANI.FR)  >>".Z) 

CALL  FSWAP( "SCAN  I 

SV") 

59 


ACCEPT  "DO  YOU  WANT  RAW  DATA  STORED?  (1=YES,  0  =  NO>  ",IRAV 
I DT ( 9 )  =  I  RAW 
IDT( 10)  =  I  PROF 
C 

C  BEGIN  TAKINC  DATA 
C 

N  =  0 
57  CALL  DI0IKISUB1  ,IRAY1  ,IRAY2,MIM,IER) 
C 

C  DECODE 
C 

AVE1=0. 
AVE2=0. 
60  DO  90  1=1, NUM 

DO  65  J=l , 16 

1B1N1 (J)=IAND(IRAY1 ( I ) ,2**<J-1 ) ) 
IBIN2(J)=IAND( IRAY2C 1  )  ,2**<J-1  ) ) 
IF( IBIN1 (J) .NE.O)  1 D 1 N 1 ( J ) = 1 
IF(IBIN2(J) .NE.O)  1B1N2(J)=1 
65     CONTINUE 
IBN1=0 
IBN2=0 

DO  70  J=0,9 

IBN1=IBN1+IBIN1 (J+l )*2**J 
IBN2=IBN2+IBIN2(J+1 )*2**J 
70     CONTINUE 
I  EXP  1=0 
1EXP2=0 

DO  75  J=10,13 

IEXP1=IEXP1+IBIN1 (J+l )*2**<J-10) 
IEXP2=1EXJ>2+IBIN2(J+I  )*2**(J-10) 
75  CONTINUE 

C 
C 

VEL1 (I )=C0NI*(2.**(5-IEXP1 )*1 000 . /FLOAT < IBN1 )-40.+FREGI  ) 
VEL2 ( I ) =C0N2* ( 2 . ** ( 5- 1 EXP2 ) *  1 000 . /FLOAT ( I BN2 ) -40 . +FREQ2 ) 
AVE1=AVE1+VEL1 (I ) 
AVE2=AVE2+VEL2(I ) 
C 

90  CONTINUE 
C 

AVE 1 = AVE 1 /NUM 
AVE2=AVE2/NUM 
C 

C  CHECK  VELOCITIES  ARE  IN  LIMITS 
C 

NN=N+1 

WR1TE(10,92)NN.AVE2,AVE1 
92  FORMATC IX, "LOCATION:  " , 12 ,4X, "U(AVE) =  ■ ,F9 . 5 ,4X, "V(AVE) =  ",F9.5) 
IF(AVE1 .CT.35. )  GO  TO  95 
IF(AVE2.GT. 140. )  GO  TO  95 
CO  TO  90 
95  WRITE  FREE(IO)  "BAD  DATA:  LOCATION  ",NN 

ACCEPT"RETRY  OR  SKIP  TO  NEXT  LOCATION?  (1=RETRY,  0=SKIP):  " , ITRY 
IF(ITRY.EO.0)GO  TO  90 
CALL  DIORM 
GO  TO  57 
C 

98  IFURAW.NE.  1  )  GO  TO  100 
C 

C  STORE  RAW  DATA 
C 

CALL  FOP EN (5, "COUNTER" ) 
READ  BINARY(5)  NJ 
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CALL  FCL0SE(5) 

NJ=NJ+1 

CALL  FOPEN( 5, "COUNTER") 

WRITE  BINARY(5)  NJ 

CALL  FCLOSE(5) 

NFILE=NJ 

IF(N.EQ.O)  JFILE=NFILE 

1R=MOD(NFILE,10) 

IS  =  NITLE/1O0 

1L=(NFILE-IS*10O)/10 

IFILEC 1 )='RA' 

IFILE(2)=,WD' 

IFILE(3)=(IS+60K)*4O0K+< IL+60K) 

IFILE(4)=(IR+60K)*400K 
C 

CALL  FOPEN(5,IFILE) 

WRITE  BINARY(5)  ( IDT( I ) , I = 1 , 10) 

WRITE  BINARYC5)  (VAR( 1 ) , I = 1 , 10 ) ,NUM, Y(NN) 

WRITE  BINARY(5)  ( IRAY1  (  I  )  ,  I  =  1  ,NUM) 

WRITE  BINARY(5)  ( IRAY2( I ) , I = 1 ,NUM) 

CALL  FCL0SE(5) 
C 

100  WRITE  (10,99) 

99  format ("  ","***  begin  rms  calculations  >>",z) 
urms=o. 

VRMS=0. 

UVRI1S  =  0. 

DO    110    1=1, NUM 

URMS=URMS+(  (VEL2(  I  )  )**2-(  AVE2  )**2) 

VRMS=VRMS+( (VEL1 ( I )**2-(AVEl )**2) ) 

UVRIiS  =  UVRM3+(VELl  (  I  )*VEL2(  I  )  )-( AVE1  *AVE2  ) 
110    CONTINUE 

URI1S=  (UllMS/FLOAT(NUM)  ) 

VKMS= (VRMS/FLOAT(NUM) ) 

UVRI1S=(UVRMS/FL0AT(NUM)  ) 

WRITE    (10, 1 1 1 ) 
1  1  1    FORMAT  ("+",">>    END   RMS   CALCULATI ONS "  ) 

N  =  N+I 

V( 1 ,N)=AVE2 

V(2,N)=AVE1 

V(3,N)=URMS 

V(4,N)=VRI1S 

V < 5 , N ) =U VRHS 
lF(N.GE.NHAX)    GO  TO    120 
C 

C   MOVE  TRAVERSE 
C 

CALL   FSWAP ( "TRELAY . SV " ) 

ACCEPT "TRAVERSE  MOVED... ANY  NUMBER  TO   CONTINUE.       ",IAA 
C 

CALL   DI0RI1 

GO  TO   57 
C 

C   OBTA I N   TEMP ,    STAG   PRES . ,    MACH   # ,    ETC 
C 

120    IF( IRUN.EQ. 1C00)    GOTO    127 
CALL   FSWAP ( "TRELAY .  SV "  ) 
WRITE(10,126) 

126  FORMAT ( "TRAVERSE   RESET   TO    INITIAL   POSITION") 

127  CALL   F0PEIU5, "SCAN I" ) 

READ    BINAIIY(5)    VAIK  4  )  ,  VAU(  5  )  ,  VAR(  6 )  ,  VAR(  7  )  ,  VAR(  8) 

CALL    FCL0SE(5) 
C 
C   WRITE    DATA    INTO   PERMANENT    DATA    FILE 
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CALL  F0PEN(5, "PERHCOUNT" ) 

READ  BINARY(5)  MF 

CALL  FCL0SE(5) 

WRITE  KREE(IO)  "LAST  DATA  FILE  USED:  "  ,MF 

WRITE  FREE(IO)  "HUN  NUMBER:  ",IRUN 

ACCEPT  "'DATA*  FILE  NUMBER  (NEC  '  TO  SKIP):  ".MFILE 

IF(MFILE.LT.O)  CO  TO  599 

CALL    FOPEN(5,  "  PEHMCOUN  T"  ) 

(rtUTE   BlNAnY(5)    MFILE 

CALL    FCLOSE(D) 

IR=MOD(MFILE, 10) 

1L=HFILE/10 

IFILE(1)='DA' 

IFILE(2)='TA' 

I F I LE ( 3 ) = ( I L+60K ) *400K+ ( I R+60K ) 

IFILE(4)=0 

CALL  F0PEN(5,IFILE) 

WRITE  B1NARY<5)  (  I  DT(  I  )  ,  I  =  1  ,  10  ) 

WRITE  BIHARYC5)  ( VA1K  1  )  ,  I  =  1  ,  1 0 ) ,NUM, NMAX, JFI LE ,NFI LE 

WRITE  BINARY(5)  ( ( V( I  ,  J  )  ,  J=  1  , NMAX) , I  =  1 , 5) 

WRITE  BINARY(5)  ( Y( I  )  ,  I  =  I  , NHAX) 

CALL  FCL0SE(5) 

599  WRITE  FREE(IO)  "AMBIENT  TEMP: 
WRITE  FREE(IO)  "ATMOSPHERIC  PRESSURE: 
WRITE  FREE(IO)  "STAGNATION  PRESSURE: 
WRITE  FREE(IO)  "FREE  STREAM  MACH  NUMBER: 
WRITE  FREE  (10)  "FUEL  STREAM  VELOCITY: 
DO  650  1=1 ,NMAX 
Y( I ) =Y( I )/6. 

650  continue 

accept "do  you  desire  a  curve?  (1=yes.  0=n0) 
if(icurve.eo.o)  co  to  973 
write  free( 10) 

600  type" what  profile  do  you  desire?" 
type"     (1)  u ( ave )/u( infinity)" 
type"     (2)  v(ave)/u( infinity) " 
type"     (3)  u(rms)/(u(infinity)**2)" 
type"     (4)  v(1u1s)/(u(ihfinity)**2)" 

accept"     (5)  uv(1uis).'(u(  infinity)**2)  ",isel 
bicj^o. 

DO  7i;0  1  =  1,  NMAX 

IF(1SEL.LE.2)  X( I )=V( ISEL, I )/VAR(0) 

IF( ISEL.GT.2)  X( I )=V( ISEL, I )/ ( VAR(8)**2 ) 

IF(X( I ) .CT.BICX)  BIGX=X(1) 
700  CONTINUE 

CALL  FOPEN (5, "PROFILE" ) 

WRITE  BINARY(5)  NMAX, I  PROF , ISEL , BIGX 

I7RITE  BINARY(5)  (X(  I  )  ,  I  =  I  ,  NM\X) 

WRITE  BINAUY(5)  ( Y( I ) , I = 1 , NMAX) 

CALL  FCL0SE(5) 

CALL  FSWAP ( " PLOTP . SV "  ) 

ACCEPT"ANY  *   TO  CONTINUE  " ,ANY 

TYPE  "DO  YOU  DESIRE  ANOTHER  GRAPH?  " 

ACCEPT"  (1=YES,  0=N0):   ",IGR 

1F( IGR.EQ. 1 )C0  TO  600 
975  WRITE  FREE( 10) "ENSURE  PRINTER  IS  ON" 

WR I TE ( 1 2 , 980 ) I  DATE , 1 T I  ME , I  V AL 
900  FORMAT(10X,"DATE:  " ,2 ( 12 ,"/"), 14 . 5X, "TIME:  " ,2( 12 , " : " ) , 12 , IX, A2) 

WHITE  FREE( 12) 

WRITE  FREE(12)"HUN  NUMBER:  ",IRUN 

WRITE  FREE(12) 

I F ( I  PROF . EQ . 1 )  WR I TE  FREE (12)" WAKE  PROF  1 LE " 


,VAR(4) ." 

DEGF 

,VAR(5) , " 

PS  I" 

,VAR(6) , " 

PS  I  " 

,VAR(7) 

" ,VAR(8) , 

"    M/S 

) :     " , I  CURVE 
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IF(  I  PROF. EQ. 2)  WRITE  FREE(  12  )  "HOUND  ARY  LAYER  PROFILE" 

WRITE  FREE( 12)"    X-LOCATION:  ",VAH(3) 

WRITE  FREEU2) 

WRITE  FREE( 12) "NUMBER  OF  Y-LOCATIONS  TAKEN:  " , NMAX 

WRITE  FREEU2) 

WRITE(12,981 ) 
981  F0RNAT(2X, "LOCATION" , 10X, "U(AVE) " , 10X, "V(AVE) " , 10X, "U(RMS) " , 
81  IX,  "V(RMS)"  ,1  1X,"UV(RMS)") 

DO  935  1=1, NMAX 

WRITE(12,984)Y(I ),V(1,I),V(2,I),V(3,I),V(4,I),V(5,I) 
984  FORMAT ( 2X , F8 . 6 , 9X , F9 . 4 , 8X , F9 . 4 , 7X , F 1 0 . 4 , 7X , F 1 0 . 4 , 7X . F 1 0 . 4 ) 
983  CONTINUE 

WRITE  FREE(12) 

WRITE  FREE( 12) "NUMBER  OF  SAMPLES  PER  Y-LOCATION:  ",NUM 

WRITE  FREE(  12)  "MIX  FHEQ  *  1  =  ".FRF.QI,"   MIX  FREQ  #2=  ".FREQ2 

WRITE  FREE(12)"VEL  CONST  #1=  ",CONl,"   VEL  CONST  #2=  "  ,C0N2 

WRITE  FREEU2) 

IFCIRAW.EQ.0)  WRITE  FREE(12)"RAW  DATA  HAS  NOT  BEEN  STORED" 

IF( IRAW.EQ. 1 )  WRITE  FREE(12)"RAW  DATA  IS  STORED  IN  'RAWD'  FILES 
8JFILE,"  THRU  " .NFILE 

WRITE  FREE( 12) "PERMANENT  DATA  FILE:  DATAH,MFILE 

WRITE  FREE(12) 

WRITE  FREE( 12) "AMBIENT  TEMP=  " ,VAR(4) , "  DECF" 

WRITE  FHEE( 12) "ATMOSPHERIC  PRESSURE=  ",VAR(5),U  PS  I " 

WRITE  FREE( 12) "STAGNATION  PRESSURE=   ",VAR(6),"  PSI " 

WRITE  FREE(  12)  "FREE  STREAM  MACII  #=    ",VAR(7) 

WRITE  FREE(12)"FREE  STREAM  VELOCITY=  ",VAR(8),U  M/S" 
999  STOP 

END 
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G  *************************** 
C  *************************** 

C  **  ** 

C  **  "SCANI .FR"        ** 

G  **  ** 

G  **  MIKE  FOREMAN       ** 

G  **  ** 

G  *************************** 
G  *************************** 

G 

G  "SCANI. FR"  TAKES  THE  DATA  FROM  THE  6 

G  SCAN  I  VALVES  AND  WRITES  TUIS  DATA  (IN  COUNTS) 

G  ON  THE  SCREEN  AND  INTO  A  DATA  FILE. 

G 

DIMENSION  IRAY(1200),  Y<6,0:23) 
G 

C  THE  SCAN [VALVES  ARE  CONNECTED  TO  CHANNELS  1  THROUGH  6: 
G 

NSTART  =  1 

LIMIT   =  6 

NUM  =  LIMIT-NSTART+1 
G 

G  ZEROIZE  THE  MATRIX 
G 

DO  10  I  =  1,6 

DO  10  J  =  0,23 

Y(I  ,J)  =  0 
10  CONTINUE 
G 

C  SCAN  I  VALVE  RELAY  INFORMATION 

C  I  DEVICE  =  DEVICE  NUMBER  (25K  FOR  SCAN  I  VALVE) 
G  RELAY  NUMBERS:  HOME  =  RELAY  *3 
G  STEP  =  RELAY  -*2 

C  TO  CONVERT  RELAY  NUMBER  TO  BINARY  WORD  USE: 
C      N=2**(16-N) 
C 

I  DEVICE  =  25K 

Nl=2        ;STEP  =  RELAY'2 

1STEP  =  2**( 16-N1 ) 

N2  =  3        ;  HOI  IE  =  RELAY-* 3 

I  HOME  =  2**<16-N2) 
G 

C  RESET  ALL  BITS  TO  ZERO 
G 

CALL  DOUTW(IDEVICE.O) 
G 

C  SETUP  THE  A/D  DEVICE  Willi  THE  STARTING  AND  ENDING  CHANNELS 
G 

CALL  FADSTC  NSTART,  LI  MIT,  I  Ell) 
C 

G  RESET  THE  SCAN  I  VALVE  TO  THE  HOME  POSITION 
C  (ENERCIZE  THE  RESPECTIVE  RELAY,  THEN  DEENERGIZE  IT) 
G 

CALL  DOUTWdDEVICE.inOME) 

CALL  FDELY(5) 

CALL  DOUTW ( I DEV I CE , 0 ) 
G 

G  TAKE  200  SAMPLES  AT  EACH  PORT 
C 

DO  30  L=0,23 

DO  20  N=l ,200 

CALL  FADCK  IRAY(  (N-l  )*NUM+1  )  .NUM.IER) 
20  CONTINUE 
G 
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C  STEr  THE  SCAN  I  VALVES  TO  THE  NEXT  PORT 
C 

CALL  DOUTV CI  DEVICE, I  STEP) 

CALL  FDELYC5) 

CALL  DUUTW(IDEVICE.O) 
C 

c  compute  the  total  counts  for  each  port  (of  each  scan i valve ) 

C  AND  STORE  EACH  TOTAL  IN  THE  ARRAY,  Y. 
C 

DO  30  n=NSTART,LlMIT 

H1=M-HSTART+1 

DO  30  N=l ,200 

Y(M1 ,L)=Y(M1 ,L)+IRAY( (N-l )*NUM+ CM-NSTART+1 ) ) 
30  CONTINUE 
C 

C  RETURN  THE  SCAN  I VALVES  TO  THE  HOME  POSITION 
C 

CALL  DOUTV  (  I  DEV  I  CE  .  I  DONE  ) 

CALL  FDELYC5) 

CALL  DOUTVC IDEVICE.O) 
C 

C  DISPLAY  THE  AVERAGE  SCANIVALVE  VALUES  (IN  COUNTS)  ON  THE  SCREEN 
C  NOTE:  THE  AVERAGE  IS  COMPUTED  INSIDE  THE  WRITE  STATEMENT  ITSELF 
C 
C 

WRITE  FREEC10)"         BEGIN  TEMP.SV..." 

CALL    FSVAP  (  "TEilP  .  SV    ) 

WRITE   FREECIO)"  END  TEIIP.SV    >>>    BEGIN   CALPRES  .  SV  .  .  .  " 

CALL  FSVAP ( " CALPRES. SV" ) 

I/RITE  FREE(IO)"         END  CALPRES.  SV  >>>  BEGIN  NPARO.SV... 

CALL  FSWAP C ■ NPARO . SV " ) 

WRITE  FHLE(IO)"         END  NPARO.SV  >>>  CONTINUE  SCANI.SV.. 

CALL  F0PENC5, "TEMP" ) 

READ  BINARYC3)  TEMP 

CALL  FCL0SEC5) 

CALL  FOPEHC5, "CALPRES" ) 

READ  B1NARVC3)  CAL 

CALL  rCL0SE(5) 

C  \LL  I  OPEN ( 3 , " ATMPRES " ) 

HEAD  BINARY(S)  ATM 

CALL  1CL0SL(5) 

ATM=  ATM;  3377/6393        -.CONVERT  TO  PS  I 

'RITE  IREEC 10) "CAL" , CAL 

TEHPR=TEHP+459.6         ; CONVERT  TO  RANKINE 

T0TFRES=(Y(1  ,2)-Y(  1  ,0)  )*CAl./(  Yl  1  .  1  )-Y(  1  .0)  )+ATM 

STATPRES=CYC4,2)-YC4.0) ) .CAL/(Y(4. I )-Y(4,0) ) +ATM 

FMAC11  =  SORT  CABS  (  3  .  *  C  C  ( STATFRES/TOTPRES  )  **  (  -2  .  /7  .  )  )  - 1  .  )  )  ) 

TRY1 =  5.*C  C (STATPRES/T0TPRES)**C-2./7. ) )-l . ) 

TRY2=ABSCTRY1 ) 

IFCTilVl  .EQ.1RY2)    GO   TO    40 

WRITE    FREEC  10)  "STATPRES    IS    GREATER  THAN   TOTPRES    !•" 
40   Tl=TEMPR/(  1+CFMACH**2.  )/3.  ) 

A=SORT( 1 .  4*1716.*T1 )*.304G 

FVEL=A*FMACH 

CALL  F0PENC5, "SCAN I " ) 

VRITE  BINARYC5)  TEMP.  ATM,  TOTPRES.  FMACH ,  FVEL 
CALL  FCL0SEC5) 
C  REMOVE  THE  DEVICES 
C 

CALL  DIORM 
CALL  FADRM 


C 


STOP 
E^U 


65 


Q     *************************** 

c   *************************** 
c  **  ** 

G    **  "TEMP.FR"  ** 

C   **  MIKE   FOREMAN  ** 

C   **  ** 

Q     *:,':************************* 
Q    *************************** 

G 

G  "TEMP.FR"  COMPUTES  THE  CURRENT  ATMOSPHERIC 

G  TEMPERATURE  AT  THE  INLET  TO  THE  INDRAFT 

G  TUNNEL.   THIS  VALUE  IS  OUTPUT  TO  THE  DISPLAY 

G  SCREEN  AND  ALSO  WRITTEN  INTO  A  FILE  CALLED 

G  "TEMP". 

G 

DIMENSION  IRAYC200) 
G 

G  THE  TEMPERATURE  PRODE  INTERFACES  WITH  CHANNEL  8 
G 

NSTA11T    =    8 
LIMIT    =    8 

NUM  =  LIMIT-NSTART+1 
G 

G  SET  UP  THE  DEVICE 
G 

CALL  FADST(NST\RT,LIMIT, IER) 
G 

G  TAKE  200  SAMPLES 
G 

DO  10  1=1 ,200 

CALL  FADCTC IRAY( I ) ,NUM, IER) 
10  CONTINUE 
G 

G  COMPUTE  THE  TOTAL  NUMBER  OF  COUNTS  FOR  THE  200  SAMPLES 
G 

TEMP=0. 
DO  20  1  =  1  .200 
TEMP=TEMP+1RAY(I  ) 
20  CONTINUE 
C 

G  CALCULATE  THE  AVERAGE  OF  THE  200  SAMPLES 
G 

TEMP=TEMP/200. 
G 

G  THE  TEMPERATURE  SIGNAL  ENTERS  THE  AMPLIFIER  IN  VOLTS. 
G  THE  AMPLIFIER  GIVES  THE  SIGNAL  A  GODB  GAIN  (FACTOR 

G  OF  1000).  THEREFORE.  THE  SIGNAL  GOES  TO  THE  A/D  CONVERTER  WITH  DIMENSIONS 
G  OF  MILLIVOLTS.   AFTER  CONVLHSION.  THE  FACTOR  163B.3  COUNTS/MILLIVOLT 
C  IS  NECESSARY  TO  CONVERT  THE  COMPUTER  UNITS  (COUNTS)  BACK  TO  MILLIVOLTS: 
G 

TEMP=TEMP/1638.3  ; CONVERT  COUNTS  BACK  TO  MILLIVOLTS 

G 

G  THE  "TEMP"  HAS  NOW  BEEN  CONVERTED  TO  MILLIVOLTS 
G  HOWEVER.  THE  TEMPERATURE  EQUATION  REQUIRES  THAT 
C  THE  "TEMP"  BE  IN  MICHOVOLTS: 
G 

TEMP=TEMP*1000. 
G 
G 

G  BECKMAN  COEFFICIENTS 
G 

C0=3. 1922707E+01 
C1=3.5497620E-02 
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C2=-2.4786653E-07 
C 

C  CALCULATE  TEMPERATURE 
C 

TEMP=C0+C1*TEMP+C2*TEMP**2 
C 

C  WRITE  TEMP  TO  SCREEN  AND  FILE  "TEMP" 
C 
C      WRITE  FREE  (10) 

WRITE  FREE  (10)  "TEMPERATURE  (DEGF)=  " ,TEMP 

CALL  FOPEN  (5, "TEMP") 

WRITE  BINARY  (5)  TEMP 

CALL  FCLOSE  (5) 
C 

C  DISCONNECT  THE  DEVICE 
C 

CALL  FADRM 
C 

STOP 

END 
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c  *************************** 
G  *************************** 

C  **  ** 

C  **      " CALPRES. FR"       ** 

C  **  ** 

C  **      MIKE  FOREMAN       ** 

G  **  ** 

G  *************************** 

C  *************************** 

G 

C  "CALPRES.FR"   THIS  PROGRAM  COMPUTES  THE 

G  CALIBRATED  PRESSURE  TO  THE  SCAN  I  VALVES. 

G 

DIMENSION  IRAY(200) 
G 

G  THE  CALIBRATION  PRESSURE  IS  CONNECTED  TO  CHANNEL  7. 
C 

NSTA11T    =    7 
LIMIT    =    7 

NUM  =  LIMIT-NSTART+1 
C 

C  SET  UP  THE  A/D  DEVICE: 
C 

CALL  FADST(NSTART,LIMIT, IER) 
C 

C  TAKE  200  SAMPLES  AT  CHANNEL  7 
C 

DO  10  1=1, 200 

CALL  FADCT(  IRAYd  )  .NUM.  IER) 
10  CONTINUE 
G 

C  COMPUTE  THE  AVERAGE  OF  THE  200  SAMPLES 
C 

Y=0. 

DO  20  1=1 ,200 
Y=Y+FLOAT(  IRAYd  )) 
20  CONTINUE 
Y=Y/200. 
C 

G  CALIBRATION  PRESSURE  IN  INCHES  OF  HG  =  VOLTS* 10 
G 

CAL=(Y/1630.3)*10. 
C 

C  CALIBRATION  PRESSURE  IN  PS  I  =  INCHES  OF  HG  *  (3377/6895) 
G 

CAL=CAL* ( 3377 . /6895 . ) 
C      WRITE  FREE(IO)  CAL 
C 

C  WRITE  THE  CALIBRATION  PRESSURE  TO  THE  SCREEN  AND  TO 
C  A  TEMPORARY  FILE  CALLED  "CALPRES". 
C 

C      WRITE  FREE(IO) 
C      WRITE  (10,25)  CAL 

25  FORMVT  ("CALIBRATION  PRESSURE:  H,F7.4,"  PSI") 
G 

CALL  FOPEN  ( 5 , "CALPRES" ) 
WRITE  BINARY  (5)  CAL 
CALL  FCLOSE  (5) 
G 

C  REMOVE  THE  DEVICE 
C 

CALL  FADRM 

STOP 

END 
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C  "PLOTP.FR" 
G  PROFILES. 
C 


IS  CALLED  BY  "SPOILEH.FR"  TO  PLOT  SELECTED  VELOCITY 


DIMENSION 
DIMENSION 
DIMENSION 


X(25),  Y(25) 
IST1 (3) , IST2(8) 
IST5(8) , IST6(9) 


IST3C8) , 
1ST7C21 ) 


IST4(0) 
, I ST8 (31) 


CALL 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 


CHRSIZ(4) 
ISTl/89,47 
IST2/85.40, 
IST3/86.40 


67/ 

65,86,69,41 ,47,85/ 

65,86,69,41 ,47,85/ 

I ST4/85 , 40 , 82 , 77 , 83 , 4 1 , 47 , 05/ 

I ST5/86 , 40 , 82 , 77 , 83 , 4 1 . 47 , 85/ 

I ST6/85 , 86 , 40 , 82 , 77 , 83 , 4 1 , 47 , 85/ 

1 ST7/87 , 65 , 75 , 69 , 32 . 86 , 69 , 76 , 79 , 67 , 73 , 84 , 89 , 32 , 
*80. 82, 79, 70, 73, 76. 69/ 

DATA  I ST8/66 . 79 , 85 , 78 , 68 , 65 , 82 , 89 , 32 , 76 . 65 , 89 , 69 , 82 , 32 , 
*86 . 69 , 76 , 79 , 67 , 73 , 84 , 89 , 32 , 80 , 82 , 79 , 70 , 73 , 76 , 69/ 


CALL 

READ 
READ 
READ 
CALL 


F0PEH(5, "PROFILE" ) 
BlNARY(o)  NMAX,  I  PROF 

(X(I  ) 

(Yd  ) 


ISEL, 
1  =  1.  WlAX ) 
1  =  1  .NIIAX) 


BICX 


B1NARY(5) 

BINARY(5) 

FCL0GE(5) 
IF(BICX.LT. .5)  BIGX= .5 
WRITE  FRLE( 10)Y( 1 ) ,Y(NMAX) ,BICX 
ACCEPT "ANY  NUIR3LR  TO  CONTINUE: 


IRAND 


CALL  INITT(960) 

CALL  BIN  ITT 

CALL   CdilSlZ(3) 

CALL   TEIIIK3,  1023) 

CALL    SLIHXC  1  GO,  900) 

CALL   SLIIIY(  150,650) 

CALL   XFRIK3) 

CALL   YIIUK5) 

CALL   XTICS( 10) 

CALL   YriCS( 10) 

CALL    DLirrX(G.O,BICX) 

CALL    DLIMY(Y(1 ) ,Y(NMAX)  ) 

C'.LL    HPTS(NfL4X) 

CALL   CliECiaX.Y) 

CALL   DSPLAY(X,Y) 

CALL    C11HSIZ(3) 

CALL   HOVAGS( 10,525) 

CALL  VLABEL  (3, IST1 ) 

CALL  M0VABS(300,725) 

IF( IPROF.EQ.l )  CALL  HLABEL(21 

IF( I  PROF. EQ. 2)  CALL  HLABEL(31 

CALL  M0VABS(420,75) 

IFdSEL.EQ.  1  )  CALL  HLABEL( 8,1  ST2) 

CALL 

CALL 

CALL 

CALL 


IST7) 
IST8) 


IFdSEL.EQ 

I F ( I SEL . EQ 

IFdSEL.EQ 

IF( ISEL.EQ 

CALL 

CALL 

STOP 

END 


2) 

3) 

4) 

5) 
FINIT1(0,750) 
CURS  I Z( 4) 


HLABEL(8,IST3) 
HLABEL(8,IST4) 
HLABKL(8, IST5) 
1ILADEL(9,  IST6) 
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c  ******************************* 
c  ******************************* 

c  **      histocram.fr      ** 

C  **         MIKE  FOREMAN        ** 

C  **  ** 

(2   ******************************* 
C  ************** ***************** 

c 

G  histogram.fr  records  data  from  the  ldv 

c  counter  and  plots  a  histogram  of  this  data. 

G 

C  REQUIRED  SUBROUTINES:  DI022  feFLIBG 

G 

DIMENSION  1RAYC 1000) , IDIN( 16) .FREQ(IOOO) 
DIMENSION  X(  101  )  ,Y(  101  )  ,VKL.(  1000) 

TYPE  "INPUT  SUBDEVICE  CODE  TWO  DIGITS  ONLY:" 

TYPE  "  CHANNEL  1  =46" 

TYPE  "  CHANNEL  2  =  47" 

READU1.5)  ISUB 
5  FORMAT (02) 

ACGEPT  "VELOCITY  CONSTANT?  " ,CON 
AGGE1T  "MIXING  FREQUENCY?  ".FREQM 
10  ACCEPT  "NUMBER  OF  SAMPLES  (0  TO  STOP)  ?",NUM 
IF  (NUM.CT. 10OO.OR.NUM.LT.O)  GOTO  10 
IF  (NUM.EQ.O)  GO  TO  999 
CALL  DIOII  (  ISUB,  IRAY.NUM,  I  Ell) 
C 

DO  100  J=l ,NUM 
DO  50  1=1,16 

1BIN( I )=IAND( IRAY(J) ,2**(I-1 ) ) 
50     IF( IBINC I ) .NE.O)  IBIN(I)=1 
IBN  =  0 

DO  60  1=0,9 
IBN=IBN+IBIN( 1  +  1  )*2**I 
60     CONTINUE 
C 

IEXP=0 

DO  70  1=10,13 

I  EXP  = I EXP+ 1 B I N ( I  + 1 ) *2** (1-10) 
70     CONTINUE 
G 

FREQ<  J ) =32000 . /FLOAT ( I BN*2** I  EXP ) 
VEL ( J ) =CON* ( FREQ ( J ) -40 . +FUEQM ) 
100  CONTINUE 
G 

AVER=0. 
AVE=0. 

DO  110  1=1, NUM 
AVER=AVER+FREQ< I ) 
AVE=AVE+VEL( I ) 
1 10  CONTINUE 

AVER=AVER/NUM 

AVE=AVE/NUM 

WRITE  FREE( 10)  "AVERAGE  DOPPLER  FREQUENCY=  ", AVER, "MHZ" 

WRITE  FREE(IO)  "AVERAGE  VEL0CITY=  ",AVE,"  M/S" 

ACGEPT  "DO  YOU  WISH  A  HISTOGRAM?  (YES=1,  N0=0)",AA 

IF(AA.EQ.O)  GO  TO  999 

MM=IFIX(AVE) 

JJ=MM-5 

IF(.J.I.LT.O)  JJ  =  0 

KK  =  J.I+10 

LL=(JJ*10)-1 
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IW=0 

Z=0. 

DO   200    1=1,101 

K=I+LL 

Xd )=(FLOAT(K))/10. 

Y(l  )=0. 

DO    150   J=l ,NUM 

I I=IFIX(10.*VEL(J)) 

IF<  II  .EQ.K)    Yd  )=Y(I  )  +  l 

150   continue 

WRITE  FREE  (10)  X(I),Yd) 
IW=IW+IFIX(Y<  I  )  ) 
IF(Y( I ) .GT.Z)  Z=Y( I ) 
200  CONTINUE 

Z=FL0AT( ( IFIX(Z/100. )*100)+100) 

IFCZ.GT. 1000. )  Z=1000. 

AB  =  FL()AT(JJ) 

AC = FLOAT (KK) 

VRITE  FREE  (10)  AB.AC.Z 

CALL  FOPEN  <5,"HTST0") 
WRITE  BINARY  (5)  AB,AC,Z 
WHITE  BINARY  (5)  (X(I),  1=1,101) 
WRITE  BINARY  (5)  (Yd),  1  =  1,101) 
CALL  FCLOSE  (5) 
CALL  FSWAP  (  "HISTOGRAPH.SV") 
WRITE  FREE( 10) 
WRITE  FREE( 10) 
WRITE  FREE( 10) 
WRITE  FREE( 10) 

WRITE  FREE  (10)  "  TOTAL  SAMPLES  PLOTTED:  " , IV 

CALL  D10RM 
999  STOP 
END 
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C  "H 
C  A 
G 


ISTOCRAPH.FR"   USED  BY  HISTOCRAM.FR  TO  PLOT 
HISTOGRAM  OF  LASER  DATA. 


DIMENSION  X(I0l ) ,Y( 101 ) 

DIMENSION  ISTK7),  IST2(21),  IST3U5) 


DATA  I  ST 1 /03 , 65 , 77 , 00 , 76 , 69 . 83/ 
DATA  I ST2/72 , 73 , 03 , 04 . 79 , 7 1 . 82 . 65 . 77 , 
*32 . 79 , 70 . 32 , 76 , 68 . 86 , 32 , 60 , 65 , 04 , 65/ 
DATA  IST3/70,82,69,81 ,85,69,70,67,89,32,40,77,72,90,41/ 

CALL  F0PEN(5,  "HTSTO"  ) 
READ  BINAKY  (5)  AB.AC.Z 
READ  BINARY  (5)  (X( I ) , I = 1 , 1 01 ) 
READ  B1NAIIY  (5)  ( Y(  I  )  ,  I  =  I  ,  1 0 1  ) 
CALL  FCLOSE  (5) 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
STOP 
END 


INITT(960) 
BINI1T 
TEHN(3,1023) 
VDAHSTC 15.6,3) 
SLIHX( 100,900) 
SL1I1Y(  150,650) 

xi  urn  2) 

YIIUK5) 
XT1CS( 10) 
Y1TCS(  10) 

xriTcs  (io) 

YMTCS(2) 
DLIMX(AB.AC) 
DLIMY(O.0,Z) 
NPTS( 101 ) 
CIILCK(X.Y) 
DSPLAY(X,Y) 
MOVABS( 10,525) 
VLABEL(7, IST1  ) 
MOVABS(350,725) 
HLABEL(21  , IST2) 
M0VABS(420,75) 
IILABEL(15,  IST3) 
M0VABS(900, 150) 
DRWARS < 900 , 650 ) 
C1IRS1Z(4) 
F1NI  IT(0,750) 
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